Bosonization approach to the mixed-valence two-channel Kondo problem 
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We present in detail the bosonization-refermionization solution of the anisotropic version of the 
two-channel Anderson model at a particular manifold in the space of parameters of the theory, where 
we establish an equivalence with a Fermi-Majorana biresonant-level model. The correspondence is 
rigorously proved by explicitly constructing the new fermionic fields and Klein factors in terms of 
the original ones and showing that the commutation properties between original and new Klein 
factors are of semionic type. We also demonstrate that the fixed points associated with the solvable 
manifold are renormalization-group stable and generic, and therefore representative of the physics of 
the original model. The simplicity of the solution found allows for the computation of the full set of 
thermodynamic quantities. In particular, we compute the entropy, occupation, and magnetization 
of the impurity as functions of temperature, and identify the different physical energy scales. In 
the absence of external fields, two energy scales appear and, as the temperature goes to zero, a 
nontrivial residual entropy indicates that the model approaches a universal line of fixed points of 
non-Fermi-liquid type. An external field, even if small, introduces a third energy scale and causes 
the quenching of the impurity entropy to zero, taking the system to a corresponding Fermi-liquid 
fixed point. 



I. INTRODUCTION 



The ideas that motivated the introduction of the two- 
channel Anderson model started with an attempt to de- 
scribe the non-Fermi-liquid physics of the UBcia com- 
pound and other U-based heavy fermionsi The model 
can be thought of as a mixed- valence model that, in cer- 
tain regimes, has the two-channel Kondo model^ as its 
low-energy theory. In that regard, the relation between 
the two models is as in the case of the single-channel 
Kondo and Anderson models; the Anderson model cap- 
tures the local-moment physics and provides a physical 
mechanism for moment formation, while at the same 
time describes also higher temperature and degenerate- 
level regimes for which mixed valence prevails and there 
are no localized moments. The study of models that 
capture such regimes is important, since there exist a 
growing number of compounds that are believed to dis- 
play mixed valence [besides the large number of U-based 
compounds, recently, a number of new Pr-based heavy 
fermions were proposed as experimental candidates for 
the realization of quadrupolar-Kondo ground states; for 
instance. La-doped PrPba (Ref. 3)]. Hence, improving 
our comprehension of the Anderson model physics is im- 
portant for the phenomenological description of all these 
compounds4i The peculiarity of U or Pr ions compared 
to the ions in other typical heavy-Fermion compounds 
is that, in a cubic crystal field, the 5/^ configuration of 
ions or the 4/^ configuration of Pr+^ are projected 
into a non-Kramers Fa doublet with a quadrupolar mo- 
ment. Fluctuations out of this state will hybridize with 
Kramers-degenerate states and provide a competition be- 
tween magnetic and quadrupolar (i.e., flavor) moments. 
A minimal model in the case of the uranium compounds, 
that takes into account spin-orbit and crystal-field ef- 
fects, leads to modeling those two states with F3 (fiavor) 



and Fg (spin) doublets that hybridize with Fg conduc- 
tion electrons and give rise to a description in terms of 
the two-channel Anderson model4 The situation is sim- 
ilar in the case of praseodymium compoundsi^i^ 

For the case of UBcis, it was originally speculated that 
the F3 doublet would constitute the lowest-energy ionic 
configuration. This opened the possibility of a quadrupo- 
lar route to two-channel Kondo physics, in which a lo- 
cal quadrupolar moment would be screened by conduc- 
tion electrons whose spin degree of freedom would pro- 
vide the two degenerate channels.- Nonlinear suscep- 
tibility measurements disfavored such scenario)^ how- 
ever, and indicated instead the possibility of a mixed- 
valence state. ^ The latter was also supported by de Haas- 
van Alphen measurements of uranium-doped samples of 
ThBei3.— The new scenario necessitates the study of the 
full two-channel Anderson model away from its local- 
moment regimes, where a low-energy Kondo- Hamiltonian 
description is not available. Subsequent theoretical stud- 
ies of the model firmly established the persistence of 
local non-Fermi-liquid physics over its whole parameter 
regime, including at mixed valence flS^Ji For a contrast, 
one could call this the mixed-valence route to two-channel 
Kondo physics. In contradistinction, the quadrupolar 
route, -practically dismissed for uranium compounds-, 
recently acquired new relevance in the context of certain 
praseodymium compounds. Measurements on PrPba^ii^ 
PrOs4Sbi2,^'*'^^ and other Pr-based skutterudite com- 
pounds are consistent with the possibility of nonmagnetic 
ionic ground states and either onsite or offsite (antiferro)- 
quadrupolar fiuctuations. For a summarized account of 
some of these experimental results, put in the context of 
Kondo physics, see, for instance, Ref. il6j. 

In parallel, the Kondo effect became an important 
subject of study in the field of mesoscopics and there 
are many proposals and experiments for realizing two- 
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channel Kondo systems using tunable setups such as 
quantum dotSfi^^— Experimental realizations based on 
the two-channel Anderson model may be more robust 
and allow for the observation of Kondo physics at higher 
temperatures, since the Kondo scale increases exponen- 
tially as the system goes into the mixed-valence regime; 
the Kondo physics should still be observable in quantities 
like the charge fluctuations encoded in the measurements 
of capacitance lineshapes<^ 

Since the original work of Kondo, the study of quan- 
tum impurity models evolved rapidly. It quickly be- 
came evident that other methods going beyond pertur- 
bation theory were required in order to access the low- 
temperature physics of these models. The most conspic- 
uous achievements in this front were the numerical renor- 
malization group developed by Wilson^^ and the Bethe 
ansatz solution obtained independently by Andrei and 
Wiegmann.^^'^^ Another important development was the 
identification by Toulouse^^''^^ of a solvable line on which 
an anisotropic version of the model is mapped into a 
fermionic resonant-level model. These non-perturbative 
techniques are required in order to study the crossover 
regimes and allow the unambiguous identification of the 
strongly coupled fixed point of the quantum impurities. 
Numerical rcnormalization group and Bethe ansatz were 
both applied successfully to Hamiltonians that incorpo- 
rate valence fluctuation physics^^ but mappings like the 
one discovered by Toulouse remained mostly restricted 
to exchange models with fixed valence (cf. Ref. Isil ) . In 
the present work, we address this missing linki^^ 

As compared to other non-perturbative techniques ap- 
plied to quantum impurity problems r^i^i^ bosonization 
-or Coulomb gas- based mappings are complementary 
and especially valuable in that they provide us with sim- 
ple alternative ways of visualizing the physicsji^ and, 
in particular, the different crossovers. In this paper, 
we examine a mapping between the anisotropic two- 
channel Anderson impurity model and a particularly sim- 
ple biresonant-level Hamiltonian. Our work generalizes 
the results of the Emery-Kivelson mapping for the two- 
channel Kondo Hamiltonian^^"— to a more involved and 
descriptive model that contains, as well, the physics of 
charge fiuctuations. After the mapping, it becomes sim- 
pler to identify the different crossover energy scales of 
the problem and to infer the existence of a line of non- 
Fermi-liquid fixed points that governs the low-energy 
physics-^^'^^ Using a renormalization group analysis, we 
establish the generic nature of our low-temperature re- 
sults and their relevance for the original two-channel An- 
derson model. We show how to calculate all the dy- 
namical and thermodynamical quantities of interest per- 
taining to the impurity over the full range of parame- 
ters and connect explicitly all the different temperature 
regimes of the system. We rederive, in a compellingly 
compact language, all the results obtained previously 
for the model using a variety of other non-perturbative 
techniques ^ lo , 1 2 ^43 ,44 obtain a number of additional 
results for the situation when external fields are present 



and the character of the infrared fixed points is modified. 

The rest of the article is organized as follows. In 
Sec. ini we introduce the anisotropic two-channel Ander- 
son model, including possible external fields acting on 
the impurity. In Sec. IIIIl we establish the mapping, for a 
set of couplings belonging to a certain manifold, onto a 
non-interacting Fermi-Majorana biresonant-level model. 
The mapping is carried out in a bosonized language and 
particular emphasis is put on the careful treatment of 
Klein factors. In Sec. IIV[ we perform a renormalization 
group analysis of the stability of the fixed points con- 
tained in the soluble manifold. In Sec.|Vl we identify the 
three crossover scales of the model and compute ther- 
modynamic quantities in the entire temperature range 
and for arbitrary values of the external fields. Finally, in 
Sec. IVIl we provide a summary of our conclusions and 
an outlook of the applications of the mapping to other 
problems involving two-channel Anderson model physics. 

II. THE ANISOTROPIC TWO-CHANNEL 
ANDERSON MODEL 

We shall consider a generalized version of the two- 
channel Anderson model to which we add terms that 
break the rotation invariance in spin and flavor spaces. 
This is akin to the standard practice in the case of the 
single-channel Kondo model of considering an exchange 
coupling constant that acquires a different value along 
the z-axis. We denote the Hamiltonian as 

H = iJhost + ^^imp + -ffhyb + -^^ficld + -ff 3 ■ (1) 

The first term describes the dynamics of the band elec- 
trons (V'lcr where a =1,1 and a = +,— correspond to 
the spin and flavor degrees of freedom, respectively) in 
the standard approximation of a linearized band disper- 
sion around the Fermi level fixed by the normal order 
prescription, 

-ffhost = dx : V'L (x) i-ivpOx) -ipaa {x) ■■ ■ (2) 

The second and third terms contain the isolated-impurity 
contribution and the band-impurity hybridization, re- 
spectively, 

^^imp = £s Actit + £f Xaa , (3a) 
(T a 

i?hyb = V^5I[^-a^a.(0)+^L(0)Ag,] . (3b) 

Here, we have used Hubbard-operator notation to de- 
scribe the impurity degrees of freedom (Xat = |a) (&| 
where a,b = a, a and the bar stands for the complex- 
conjugate representation). As compared to slave- 
operator notation, the use of Hubbard operators auto- 
matically restricts the Hilbert space of the impurity to 
the physical one.'*'* These first three terms constitute the 
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standard two-channel Anderson model. The fourth term 
describes the coupling to external fields: 



(4) 



where, for the sake of generality, we included two fields: 
one coupling to the impurity spin and the other one cou- 
pling to its flavor. Finally, the fifth term provides the 
generalized anisotropy and can be written as a sum over 
charge (c), spin(s), flavor (/), and spin- flavor (s/) sectors 
(as will be seen below, these sectors arise naturally after 
bosonizing the model). 



^3= E ^3 



(5) 



with each term taking the form of a density-density in- 
teraction between the different impurity densities with 
the corresponding ones from the band. 



(6) 



where = [x = Q), 



(■^)= E :V'L(^)TL,.'.'^aV'(x): (7) 



and 



Q(T,q'(t' — ^aa'^aa' i ^aa.a'a' — ^aa'i aa' 



' a.a."-'aa' i ^nin-ni'rr' ~ aa.' ' aa' •, 



^3 ^3 



(8a) 
(8b) 



(here is the third Pauli matrix). The impurity densi- 
ties involved are 

X, = E ' (9a) 

a 

X/ = E "^aa , (9b) 

Xc = Xsf — E ■^<^<^ ^ E ' (9''-' 

(T a 

Since the impurity Hilbert space contains four states, we 
can define only three independent density-like operators 
apart from the identity; we find it convenient to define 
the spin-flavor density to coincide with the charge one. 



III. BOSONIZATION BASED MAPPING 

Bosonization is a well known technique that ren- 
ders more accessible the study of 1-1-1 dimensional 
models.— 1^^— When combined with refermionization 
procedures, the bosonic language can be used to find non- 
trivial mappings between different fermionic models. The 
idea is that, many times, a mapping that results highly 
nonlinear in the fermionic language can be written as a 
simple canonical transformation in terms of the bosons. 



A. Bosonization 

We shall employ the standard bosonic representation 
of the fermionic fields, 
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(10) 



in which a is a regulator that plays the role of an in- 
verse bandwidth and c/iao- are bosonic fields that describe 
the particle-hole excitations around the Fermi sea. Fi- 
nally, are the so called Klein factors, responsible for 
recovering the correct anticommutation relations among 
different fermionic species and necessary for describing 
processes in which the number of fermions changes. They 
act as ladder operators in the fermionic Hilbert spaced 
and commute with the 0's. Additionally, they satisfy the 
following algebra: 

^aCT-^Qcr — Faa-F^^ — 1 , (H^,) 

Fl^F^,,, = -F^,„,Fi^ for [aa) + (aV) , (lib) 

FaaFa'a' = ^Fa'a'Facj for (ofCr) ^ {OL o' ) (llc) 

and they obey the following (anti)commutation relations 
with the impurity Hubbard operators: 

\Faa,Xaia"\ = \Faa-,Xc'a"\ — , (12a) 
{Faa^Xfji^i^ = {Fa„,Xa'cr'} = . (12b) 

In terms of the bosons, the Hamiltonian for the band 
takes the form 



(13) 



Following Emery and Kivelson, it is natural to introduce 
a rotated basis for the bosons = c, s, /, sf) 



+ a(f>s + a(f>f + aa(f)sf) 



(14) 
(15) 



where a,a — ± when entering as multiplying factors. 
The form of iJhost remains the same in the new basis: 



-ffhost, — 



'host — 



(16) 

Notice that the Klein factors disappeared from the 
bosonized version of the first term, but they will enter 
explicitly in the hybridization term: 



V 



hyb 



'2'Ka 



XfJCli Ffy 



(17) 

Last, the terms involving only impurity operators 
stay unchanged and the terms involving exchange are 
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bosonized according to the standard prescription for den- 
sities: 

= J^X^p, = -^XMM ■ (18) 



The difficulty for studying this model is contained in 
the highly nontrivial form of the hybridization witli the 
impurity. Therefore, the strategy is to look for a canon- 
ical transformation that simplifies this term. We define 
the following generic transformation U = UcUsUfUsf 
with 



(19) 



For transforming -ffhyb, we first notice that U commutes 
with the vertex operators and the Klein factors. Thus, we 
only need to compute expressions of the form U^XaaUl- 
Using 



[X, 

[Xcra,Xf] = aX, 



ctX„ 



era 1 -^c] — 2^^Q, , 



[X^^,X,] = 0, (20a) 
[X^^,X,] = 0, (20b) 
(20c) 



(note that the two identities on the right imply that iJimp 
and -fffieid are not affected by the transformation), we 
obtain: 



UgXaaUl = XcraG 



UfXcraUj- 
UcX„^Ul 

UsfX^s^Ulj. = X^^e^'-'^f^^f 



(21a) 
(21b) 
(21c) 
(21d) 



Therefore, we can transform the hybridization Hamilto- 
nian as .ffhyb = UH\yy\Jj\ with 



lyb 



X e 




\x - F p-i^h-'^io)<t>cp-i<x(^+ij)4>j 



'e-*(^-2w)^./+h.c. 



(22) 



Next, wc use the freedom of choosing the particular trans- 
formation that will simplify the form that this term takes. 
Such a choice is given by taking 7c = 1/4, 7/ = —1/2 
and 7s = 1/2, and thus decoupling the impurity from 
the corresponding sectors in the band. The transformed 
expression is 



in the first term while giving the second one fermionic 
dimensions: 



-ffhyb 



V 



[(XT+FT++X4_i^4_) 



+ {X^_F^_ + X^+F^+) e"^'f + h.c] . (24) 



The spin-flavor band and the impurity remain coupled, 
but the form of their hybridization is now much simpler. 
Before studying this term further, let us indicate what 
happens to the remaining terms in the Hamiltonian upon 
the same transformation: 

U (ifhost + Hs)U^ = i?host 

+ - y^invFj, - J,)d^MO)X, . (25) 



We define for future reference the couplings 



(26) 



which can be set to be zero by tuning the values of —^ 
TTVFj„; thus rendering the model particularly simple, as 
it will be seen below. 



B. Refermionization 

Inspired by the classic results by Toulouse and by 
Emery and Kivelson in which they map the single- and 
the two-channel Kondo models, respectively, into differ- 
ent types of resonant-level models and find particular 
values of the anisotropy for which those resonant-level 
models are non-interacting and therefore exactly solv- 
able, we find it compelling to seek a similar mapping 
for the two-channel Anderson model. However, while 
in the Kondo model the interactions between the band 
electrons and the impurity are exchange terms that in- 
volve only fermionic bilinears, that is not the case for 
the Anderson model in which the hybridization term al- 
lows for the fiow of charge between the band and the 
impurity. This is an important difference that gives rise 
to a much richer parameter regime, in the case of the 
Anderson model, that includes not only local-moment 
phases but also mixed-valence ones. Prom the techni- 
cal point of view, it becomes important to keep track 
of the charge transfer processes while bosonizing, which 
requires a proper treatment of the Klein factors. 



Hhyh = 



V 



+ (Xt_F|_ + Xi+F^+) e-'(-5-273/)<^,/ + h.c.] , (23) 

(where we have omitted the bar over the values taken 
by the a subindex). We can further simplify this term 

by choosing 7^^ = ±1/4. The two choices are equivalent 
and we choose 75/ = 1/4; in this way, we kill the vertex 



1. New fermions 

In order to find a mapping to a new fermionic model, 
we need to refermionize the version that was already sim- 
plified by the cal transformation just discussed. The first 
step is to introduce a new set of Klein factors that allow 
us to define new fermionic operators corresponding to the 
physical sectors (c, s, /, sf): 
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(27) 



The new Klein factors are defined to be ladder operators 
in the physical Hilbert space for the band and to obey 
the usual algebra given in Eqs. piap - (|12b[) . 

Inspecting 77hyb, we are led to define the new impurity 
operators: 



d- 
f 



X 



^1+^1+) 



(28a) 
(28b) 



They allow us to rewrite the hybridization term as 
V 



hyb 



V 



/)] 



(29) 



The new form of the Hamiltonian looks compellingly sim- 
ple. However, we will find that the new operators are not 
fermions. It can be seen that they satisfy the following 
relations: 



and 







, 


(30a) 


dd'' 


= X|| - 




(30b) 


dU 






(30c) 


f 


= = 


, 


(31a) 


fP 


= 


+ -'^TT ' 


(31b) 


Pf 




+ ■ 


(31c) 



Using the completeness relation for the impurity Hilbert 
space (X||+X||+X+++X__ = 1), it is clear that both d 
and / are self -fermions {i.e., each of them independently 
obeys fermionic anticommutation relations with itself). 
We still need to verify the mixed commutation relations. 
It is important to stress that the relations written above 
can be deduced independently of any relation between 
old and new Klein factors. That is not the case for the 
relations mixing d and /. 



2. Klein factor relations 

Following the authors of Rcfs. 'Td and^A§., we have the 
freedom to make the following four identifications be- 
tween the sets of old and new Klein factors: 



fIf! = fLf, 



FsfFj 



F^F^^ 



FlFj = FT Ft, 



sf^ f 

F^F^ 



■^T+-^T- 



(32a) 
(32b) 
(32c) 
(32d) 



These relations are consistent with the physical con- 
straints given by the variations in the fermionic num- 
bers on the different sectors the Klein factors act upon, 
and the choice of signs serves to fix otherwise arbitrary 
phases. All other relations and their respective phases 
are now automatically fixed, for instance. 



FtF; = f/^Fx_ 
FsF^-F^Fl 
F}Fsf = FlF,. 



(32e) 
(32f) 
(32g) 



can be deduced using the first four relations. 

We want now to explore further the relations between 
old and new Klein factors. For example, using Eqs. (|32ap 
and p2b)l . one can derive the following expressions: 







:{Ft+,Ft^n+} = {FT+,Fj^Ft} 
Ft+Fj^{FT+,Ft} 



Ft+,F/_Fi_ 



Ft. 



= {F^+,F,f}Fj-Fsf{F^+,F}} 



(33a) 



(33b) 



Premultiplying these equations by F,/ and fJ^, respec- 
tively, and adding them up, one finds 



FsfF-^+Fj^ 



FlfFnFsf 



(34) 



Using the different relations in Eqs. (|32ap ~( P^gl ), one can 
generalize this relation to a generic relation between old 
(o) and new (n) Klein factors: 



(35) 



Interestingly, such a relation is not consistent neither 
with commutation nor with anticommutation relations 
between the old and new Klein factors. This is not unex- 
pected, since the two sets exist in different Hilbert spaces 
and the physical identification between the two is only 
through the relations among bilinears. While in Kondo- 
type models this poses no problems because only Klein 
factor bilinears enter in the Hamiltonians, this is not the 
case in Anderson-type models, '^^ for which we need to go 
beyond bilinears. Therefore, if we insist that the old and 
new Klein factors should obey a relation of the type 



F F — (tF F 



(36) 



we find that = —1 must hold. Thus, only semionic 
commutation relations are consistent between old and 
new Klein factors and there is still an arbitrariness in 
the phase that needs to be fixed. Let us postulate 



F-\+Fc — iFcF^^ 



(37) 



all other commutation relations can be deduced from this 
one and summarized as 



F F^ — —iO'^ F'f F 



(38a) 
(38b) 
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where 6'i 



a, 



eL = cr, and 



aa. 



A similar procedure can be used to verify that the new 
Klein factors and the nondiagonal impurity Hubbard op- 
erators that enter in iJhyb obey 



FX — X F — () 



(39) 



which is consistent only with standard commutation or 
anticommutation relations between the operators and, on 
physical grounds, we take them to anticommute. 

Using the semionic commutation relations between the 
two sets of Klein factors, it is easy to see that 



df = ifd , 



(40a) 
(40b) 



Analogously, with respect to the spin-flavor band, we find 



and 



[/,v./] = = 



(41a) 
(41b) 

(42) 



We say that d is a relative semion with respect to / and 
■08/; while the latter two are relative bosons. 

On the other hand, the band fermions from the charge, 
spin and flavor sectors, do not have simple commutation 
relations with the new impurity operators. Those can be 
simplified by performing a unitary transformation of the 
Klein factors belonging to that subspace. Namely, we 
define 



F, 



Ff = e^'^^fFf 



Fr 



(43a) 
(43b) 
(43c) 



The redefined operators in the decoupled bands are now 
relative fermions with d, f, and among themselves. No- 
tice also that the Hamiltonian for the decoupled sectors 
is invariant under the redefinition. What remains is to 
simplify further the commutation relations among the 
spin- flavor operators and the impurity ones. 



3. Jordan- Wigner procedure 

The procedure we will use to simplify the commutation 
relations among the impurity and spin-flavor operators is 
inspired on the Jordan- Wigner treatment of spin chains. 
We have that ipsf, d and / are all self-fermions while 
the pairs {ipsf,d) and (d, /) are relative semions and, fi- 
nally, {ipsf, /) are relative bosons. For the sake of clarity, 
we will split the Jordan- Wigner transformation into two 
parts. In the first part, we shall change all operators into 
relative bosons, and in the second part, we will perform 
a standard Jordan Wigner transformation to turn them 



into relative fermions; that way the full system will be a 
system of fermions. 
Defining 



d — d , 



(44a) 
(44b) 

Vi./ = e-""-/2i/-,/ , (44c) 



(where = d'^ d and = /^/), we see that all operators 
are still self-fermions, but now their relative statistics is 
that of relative bosons: 



[d,/] = [d,^,/] = [/,Vs/]=0. 



(45) 



A set of operators that are relative bosons but all self- 
fermions, is what is usually called core bosons. They are 
equivalent to spin- 1/2 spins and the usual Jordan- Wigner 
treatment is available to turn them into fermions. The 
second transformation is thus written as 



/ = / 



Ipsf 



^sf , 



(46a) 
(46b) 
(46c) 



(notice that we have redefined Tpsf, d and /, and that 
?!/ = n J and Ud = n^). While the first Jordan- Wigner 
string was anchored at the d site, the second one is an- 
chored at the / site; this is why it is physically more 
transparent to split the transformation into two steps. 

Finally, since the different fermionic species can be or- 
dered such that 77hyb contains only nearest-neighbor hop- 
ping terms, it can be seen that the Jordan- Wigner strings 
will not appear explicitly in the Hamiltonian. After the 
double replacement, we find 



hyb 




f){d^+d)] , (47) 



where the operators are now the new ones and we were 
left with a system where all operators are standard 
fermions (including those in the decoupled bands). 



C. Biresonant-level model 

Rewriting the remaining terms of the Hamiltonian in 
the language of the new fermionic operators and choosing 
the anisotropy that would make Aj/ = Vi', we find that 
the transformations outlined above, in the limit of zero 
fields, give a mapping to the following Fermi-Majorana 
biresonant-level model (up to an additive constant, see 
also Fig. [1]): 



biRcs 



= H'J -e dU + V2A [vl/ (0) d + dHsf (0) 

-f V2f (/t-/)(dt+d) , (48) 
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FIG. 1: Schematic representation of the mapping procedure 
between the anisotropic two-channel Anderson model and the 
biresonant-level model. The external field enters as a site 
energy for the fermion and has the effect of coupling its 
otherwise disconnected Majorana component. 

where £ = - £/, A = V"^ 12 and V = A/27ra. This 
is a purely quadratic model, on which reintroducing the 
terms with non-zero would parametrize the deviations 
from the solvable manifold (cf. Refs. '49'!50'). 

Rewriting the iJficid term after the mapping, we find 

F|gid = /t(2n/-l)-A,i(nd + n/-i) , (49a) 
= 2A,,n/nd , (49b) 

where we defined A/j = hf — lis and 2/i = hf + hs. Taking 
A/i = we find that the model is still quadratic even 
for non-zero h. Nevertheless, we will see below that the 
presence of this term changes the fixed point; we will talk 
in this case of a double-Fermi biresonant-level model. 

It is interesting to compare these two resonant-level 
models between them and with the models that cor- 
respond to the solvable points of the single- and two- 
channel Kondo models. In the latter cases, one obtains a 
single-Majorana resonant-level model and a single-Fermi 
resonant-level model, respectively. In contrast to the 
Kondo models, the Anderson case requires two fermionic 
degrees of freedom to parametrize the state of the im- 
purity, which is now more complex, since it allows for 
fluctuating states of valence. For both the two-channel 
Kondo and Anderson models, in the absence of external 



fields, "half of a fermion" coming from the impurity de- 
couples from the band and the models are naturally writ- 
ten in terms of the Majorana components of the impurity 
degrees of freedom. The Majorana component that ex- 
ists completely decoupled from the rest of the system is 
responsible for the fractional residual impurity entropy 
that we will discuss below and constitutes a signature 
of the non- Fermi-liquid properties of the fixed point. In 
both cases, the addition of an external field has the effect 
of coupling the disconnected Majorana component and 
driving the system toward a Fermi-liquid type of fixed 
point, like the one of the single-channel Kondo model. 



IV. STABILITY OF THE FIXED POINTS 

As mentioned above, taking A^; = (for x ~ c, s, /, s/, 
h) defines a solvable manifold that contains a set of fixed 
points. The important question to be addressed is how 
generic and stable are those fixed points, how are they 
parametrized and whether RG-flow trajectories starting 
outside the manifold will flow to the same fixed points. 
The last consideration is particularly important, since it 
addresses the genericalness of the fixed points and the 
feasibility of perturbative calculations in A^;. 



A. Inside the soluble manifold 

In order to address these questions, we start by com- 
puting the different Green functions. For that, we first 
write down the local action for the impurity, in which 
the extended degrees of freedom from the band were in- 
tegrated out exactly; this procedure is justified in order 
to study the interaction of the band with a local impurity 
that exists only at, say, ccimp = 0. In Nambu notation, 
we write all the sums over positive Matsubara frequencies 
(ujn) only and use the following spinor basis: 



= ( ^sf (Wn) Iplf (-t^n) d (Un) (-W„) / (w„) P {~UJn) ) 

The explicit form of the local action in the solvable manifold is 

5= ^^^'t(^„)A(^„)*(c^„) , 



n>0 



where (3 = I/ZcbT is the inverse temperature and 



A{uJn) 



/ 


—2ivF 





V2A 








° \ 







^2ivp 





-V2K 










V2A 





-iujn - e 







-V2r 







-V2A 





-iuJn + £ 


\/2r 
















—i(jJn + 2h 





V 








-^/2^ 


-V2r 





—ioJn — 2h / 



(50) 



(51) 



(52) 
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The two-point Green functions are defined as (3G{uJn) = (^(w„)^''''(w„)) with G(w„) = A ^{ujn)- Given the structure 
of the model, it is convenient to work in terms of Majorana fermions?^ For that, we rotate each spinor according to 



I 



(53) 



where i? is a block diagonal matrix in which each block, 
given by 



J_ f-i I 
V2 I 1 V ' 



(54) 



rotates a given Nambu doublet. With this conven- 
tion, the Majorana Green function matrix is GMi'^n) = 
RG{uJn)R^ and its expression in imaginary time is 



-Gl,(c.„)e*'^"(^-^') 



(55) 



Let us now discuss the situation when h = (the case 
of /i 7^ will be discussed in the Appendix). A direct 
calculation gives the following results for the long-r be- 
havior of the different diagonal Green functions: 

(V;}(r)V4(0)) ~ i , (^;^(r)V.;/(0)) ^ ^ , (56a) 



T 

(d"(r)d"(0))~i 
{f"{r)f"m ^ \ 



(d'(T)d'(O)) ~ - , (56b) 
(/'(^)/'(0)) ^signr . (56c) 



Notice that different Majorana components of a given 
fermion exhibit different asymptotics; this fact and the 
behavior of certain Green functions were already 
stressed previously in the context of the two-channel 
Kondo mode h^^i^^ From these correlators, one could read 
the following naive scaling dimensions (see also Fig. [5]): 

[€f] = K] = [/"] = 1/2, Wsf] - K] = 3/2, and [/' ] = 
0. The long time behavior of all but one of the mixed 
correlation functions can be obtained using these scaling 

dimensions; the exception being (^■ijj'^j:{T)d" (0)^ ^ l/r^. 

The scaling dimensions above are not compatible with 
a standard renormalization group scheme that proceeds 
by integrating energy shells while keeping certain coeffi- 
cients of the quadratic action constant, cf. Ref. HI. In- 
stead, we follow the ideas of Ref. [13 and introduce the 
Majorana Fermi velocities {v^i, v^n) as well as constants 
multiplying the frequency in the Berry terms for the im- 
purity {gd' , gd", gf'jgf") and let all of them vary. Specifi- 
cally, we proceed along the standard stepa^ but insisting 
that the fields should scale with the dimensions we deter- 
mined from the Green functions, this yields the following 



[w']=3/2 [d"]=l/2 
— ■< ►— 



[f']=0 



r 
— 



h=0 



[y"]=l/2 [d']=3/2 [f"]=l/2 



FIG. 2: Connectivity graph, in terms of Majorana operators, 
for the biresonant-level model. The scaling dimensions of the 
fields are seen to alternate between 1/2 and 3/2, with the 
exception of /' that stays decoupled (for h = 0) and has zero 
scaling dimension. 



RG-flow equations for the different couplings: 



dgd' 

de 


= -3gd' , 


dgd" 

de 


= -gd" , 


(57a) 


dgf 
de 


= , 


dgf" 

de 


= -9f" > 


(57b) 


dv^i 

de 


= -2w^' , 


dv^" 

de 


= , 


(57c) 


dA 

u 


= "A, 


de 

de 


= -e , 


(57d) 


dV 

m 


= -r , 


dh 

le 


= h. 
2 


(57e) 



We find that most couplings are RG-irrelevant but for 
v^" and gf that are RG- marginal. In here, we are con- 
sidering h = 0, but we find that if a field were to be 
generated it would be RG-relevant. It is important to 
remark that the the quotient e/A obeys 



deKi)=' 



(58) 



and is also marginal. This agrees with the results of 
Bethe ansatz^^'"'^ and boundary conformal field theory^ 
that find a line of non-Fermi-liquid fixed points that, in 
the microscopic theory, are parametrized by the e/A ra- 
tio and interpolate between spin and flavor two-channel 
Kondo behaviors. 



B. Outside the soluble manifold 

Collecting the terms in the rcfcrmionizcd Hamiltonian 
that takes us away from the soluble manifold, we have 

Hx = —XcPc — ^sfPsf — XsPs 

+ 2XcPcnd + 2\sfPsfnd + XfPfrid 
+ XsPsn-d + IXsPsUf 
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2XfPfndnf - 
Ah (ud + rif 



2XsPsndnf 
- i) + 2Xhnfnd 



(59) 



In order to determine the stability of the h = fixed 
points, we proceed to study the RG-relevance of the dif- 
ferent terms in Hx. The first three terms are just chem- 
ical potential terms for the conduction electrons. The 
charge and spin ones pertain only to bands that are de- 
coupled from the impurity and are marginal operators, 
since the dimension of the respective fields is 1/2, as dic- 
tated by their Berry phases. On the other hand, the spin- 
flavor one is irrelevant as can be seen by simple power 
counting. 

For the operators in the second, third and fourth lines, 
we get the following dimensions by summing the dimen- 
sions of their constituents [pc^-d] — [Pfi^d] = [Psn-d] = 3, 
[Psfnd] = 4, [psUf] = 3/2 and [pfUdUf] = [ps?^<^?^/] = 
7/2. Hence, all of them are RG-irrelevant. None of these 
terms can generate a local field term for the impurity. In 
particular, ipc, ips and ■0/ are decoupled (and the average 
of the normal-ordered densities pc, Ps and p/ is zero). 
Eventually, the operator PsfUd may give a correction to 
e, thus taking the system along the line of fixed points. 

An interesting situation appears when one starts from 
the fixed point in the absence of a magnetic field (/is = 0) 
and considers the effect of a small \h 7^ (this is realized 
when one starts with A^, = 2/i, or in other words, with 
hf ^ 0). Then one must consider carefully the effect of 
the last term in H\, with dimension [n/n^] = 5/2 and 
therefore, in principle, RG-irrelevant. This term is actu- 
ally classified as dangerously irrelevant, as it generates 
a potential field for / {i.e., an RG-relevant h term). It 
thus couples the remaining Majorana component, /', and 
changes the fixed point. 

Indeed, when adding a magnetic field the fixed point 
changes (see Appendix) . The chemical potentials for ^jjsf 
and d become marginal, and the operators with four and 
six fermions are now all irrelevant. The parameters e/A 
and h/V parametrize the new set of Fermi- liquid fixed 
points. 



V. THERMODYNAMICS WITH FINITE FIELD 

In order to extract the impurity thermodynamics when 
the model parameters are on the solvable manifold, we 
can resort to an exact calculation of the free energy. The 
impurity free energy can be straightforwardly calculated 
by means of Pauli's trick of integration over the coupling 
constants (an alternative way is to first derive an effective 
action for the impuritj*^), namely, 



(60) 



The ^ subindex indicates that the mean values should be 
computed using an action in which all couplings involv- 
ing the impurity were multiplied by the dimensionless 



parameter ^, and Qq is the impurity free energy when 
^ = 0. After computing the mean values using (wn), 
one arrives at the expression 



where 







1 



(61) 



- [AuJn (8r + Acj„) + Ah^ (e^ + 2Aw„)] 

- ujI {4h^ + 8T + e^ + 2Aw„) + . (62) 



As a function of this is a fourth order polynomial 
with roots that are parametric functions of ^: ujk 
/c = 0, . . . , 3. Using the factorized form of the polynomial 
in terms of its roots and introducing a suitable regular- 
ization, we arrive at the following expression: 



n>0 k 



(63) 



We verify that LOk{^) = and call uuk = w/c (1). Exchang- 
ing the sum in n and the integral in /x, we arrive at the 
final expression: 



k 

27r 



, 1 A/3\ 



27r 



, (64) 



where (z) = 9z InF (z) is the digamma function. Note 
that the free energy needs the A regulator, but in all the 
derived thermodynamic quantities we will be able to take 
the the limit A — > +00 and obtain finite results. 

In particular, the impurity entropy is given by S"— 6*0 = 
X]fc s (zfc), with Zk = -/3a;fe/27r and 



s{z) 



1 



-InF 



1 



hiTT. (65) 



Remarkably, the function s (z) is the same that is found 
for the entropy of the two-channel Kondo model on 
the Emery-Kivelson line.^^ This is not completely un- 
expected, since Kondo is the low-energy effective theory 
for most part of the parameter regime of the two-channel 
Anderson model. In Fig. [3l we show the impurity con- 
tribution to the entropy as a function of temperature for 
different values of the local field h. 

Let us first discuss the h — Q case shown on the 
top left panel of the figure. In this case, one of the 
roots is zero (wa = 0) and there is a residual entropy 
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FIG. 3: Impurity contribution to the entropy as a function of temperature for different values of the symmetric external field h. 
Temperature and field are both measured in units of A. Different curves correspond to different values of e/A and F — A/27r. 
The top left panel corresponds to zero field and the other ones to finite values of the field that increase hundredfold between 
panels, from h = lO'^^A until h = lO^A. 



of S'imp = fee In \/2 that signals a non- Fermi-liquid set 
of fixed points. The finite temperature physics is gov- 
erned by the other three roots. We find that, within 
the range of parameters of physical relevance, one of the 
roots (wq) is real while the other two are complex 

conjugate of each other. This makes it natural to identify 
the roots with the Kondo and Schottky energy scales in 
the following way: 

^bTk = -wo/27r , (66) 



feTs = |wi|/27r . (67) 

The Kondo scale corresponds to a jump in the entropy 
of hight fee In \/2, while the Schottky scale is associated 
with a jump that is twice as large. For small |e|, the two 
quenching steps coincide; whereas as |e| /A grows, the 
Kondo temperature decreases while the Schottky tem- 
perature increases. Evidently, over the soluble manifold, 
the dependence of all scales on the microscopic parame- 
ters is algebraic. This contrasts with the isotropic case. 
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for which the Kondo scale has exponential dependence 
on \e\ /A.'*^ This behavior is seen also in the single- and 
two-channel Kondo models, for which it is found that 
the exponential dependence of the Kondo scale is a prop- 
erty of the isotropic models only and the functional de- 
pendence crosses over to algebraic in the Toulouse and 
Emery-Kivelson limits, respectively.6S,i^ Another char- 
acteristic of the solvable limit is that the InT depen- 
dencies in the impurity susceptibilities or in the specific 
heat coefficient are changed into power laws. It is known 
that the logarithms can be recovered using perturbation 
theory in (cf. Refs. [HillO), as it reintroduces the 
leading RG-irrelevant operators that are absent from the 
soluble manifold. For the leading contribution, the three 
operators with scaling dimension 3 (i.e., pcnj,, Psnd, and 
PfUd) can be included using second-order perturbation 
theory in the couplings. All three of them mix one of the 
decoupled bands and the impurity in a way that is the 
direct generalization of what happens in the case of the 
two-channel Kondo modeli^ Such a procedure should re- 
cover, for instance, not only the results for the specific 
heat and susceptibilities but also the Wilson ratio that 
was discussed already in Ref. [4ll 

When /i 7^ 0, W3 becomes non-zero as well and grows 
with h. This introduces a third quenching step (as it 
can be seen, for instance, in the top right panel of the 
figure) in which the entropy goes from fee In y/2 to zero, 
signaling the system fiowing away from the non-Fermi- 
liquid fixed point. This is in accordance with the result of 
the previous section, in which we found /i to be a relevant 
perturbation. This process defines a third energy scale, 



(68) 



that grows with h until reaching the same value as Tk, - 
which happens at different values of the field for different 
values of |e| (taking place first for larger values, see the 
middle left panel in Fig. [3]) . As these two energy scales 
merge, the respective roots become a complex-conjugate 
pair, a single quenching step of the same height as the 
Schottky one emerges and the transition width (i.e., the 
width of the corresponding anomaly in the specific heat) 
narrows. As h increases further, the three scales first 
become degenerate and then cross each other. In the 
bottom right panel, the two scales have already crossed, 
the higher temperature step is given by Tk = and 
is insensitive to the value of |£| / A, while the lower one 
continues to correspond to the Schottky transition and 
shows the same type of dependence in |£| / A as for lower 
fields. These results share many properties with those 
for the isotropic case, but they are specific to the case of 
a symmetric filed (A/j = 0) and exhibit important differ- 
ences with the case when only hg or hf is applied to the 
system (cf. Fig. 7 in Ref. HI). 

Other thermodynamic quantities of interest are the im- 
purity charge valence Uc = dn/ds and the magnetization 
Wimp = —dQ/dh. They are explicitly given by the ex- 




FIG. 4: Impurity charge valence as a function of temperature. 
The curves are displayed for h = 0, the only effect of finite 
fields is to wash out the overshoot that happens for small val- 
ues of e/A and to shift (away from 1/2) the zero-temperature 
values of ric (see text). 



pressions: 



TO 



imp 



'T^imp.O 



^Y."^ \\ + ^hZk , (69b) 



/3- 



where the derivatives of the roots can be expressed in 
closed form using the identity 



ClE,h^k — -f 



(70) 



Plots of these two quantities are displayed in Figs. [4] and 
\5\ respectively. 

The variability of the impurity valence is an aspect 
inherent to Anderson-type models and of relevance, for 
instance, in the context of quantum dots and other meso- 
scopic systems that might allow direct measurements of 
the valence states via changes in capacitancei^ The va- 
lence starts as Ucfi = 1/2 at high temperature (T :$> Tg) 
independently of the values of e and h, and evolves 
as the temperature is lowered to attain a certain zero- 
temperature value = ric (e, h)rp^Q that characterizes 
the particular fixed point that the system reaches. The 
quenching of the valence fiuctuations coincides with the 
first quenching of the entropy and takes place at the char- 
acteristic scale Tg. Subtle aspects of the small |e| curves, 
such as the "overshoot" of the curves at intermediate 
temperatures, T < Tg, are generic and are also present 
in the exact solution of the isotropic model.— i^i^ This 
feature exists for h zero or small and disappears as h be- 
comes of the order of A. Comparisons with the isotropic 
model show that, despite the differences in the infrared 
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1.0 h ' ' ' ' 

;i/A=io"' 

0.8 - 




10-' 10-* 0-01 1 

77A 



FIG. 5: Impurity magnetization versus temperature. From 
top to bottom, the three panels show three different values of 
ft/A = 10"^10"^10"^ Curves for different values of e/A 
are shown and the values are indicated in the figure. 



physics, the isotropic and the anisotropic models share 
the same generic ultraviolet physics. "^^ This is to be ex- 
pected, since the physics at higher energies tends to be 
dominated by local fluctuations. 

The curves of finite-field impurity magnetization as a 
function of temperature are in some ways similar to those 
of the charge valence. At high temperatures, the magne- 
tization is zero, and it acquires a field-dependent finite 
value as the temperature becomes T <Th (the magneti- 
zation does not reach the maximum value of 1 due to the 
hybridization of the impurity with the conduction band). 
In this case, curves for positive and negative values of e 
are degenerate. The quantities m^-^^ = ni-imp {s, ^)t=o 



and n° label the set of finite-field fixed points of the 
model. As we mentioned above, these fixed points are 
of the local Fermi-liquid type. 



VI. CONCLUSION 

We have shown how the addition of exchange terms 
that provide anisotropy that breaks the SU (2) symme- 
try of the two-channel Anderson model in the spin and 
flavor sectors, by singling out an easy access for each of 
them, allows to identify a manifold in coupling constant 
space over which the model becomes quadratic and thus 
exactly solvable. The solvable theory is obtained in the 
form of a mapping to an anusual type of resonant-level 
model that includes two resonant levels and we thus call 
it a biresonant-level model. Similar mappings exist for 
the Kondo model in the case of one and two channels 
(and are known to be impossible for larger number of 
channels^-). But a simple mapping of this type is new 
for a model like the two-channel Anderson model that 
incorporates the physics of valence fluctuations. Equiv- 
alent mappings for the single-channel Anderson model 
are elusive (cf. Ref. Isil ). which indicates that the case 
of the two-channel model constitutes a singular exam- 
ple. Since the mapping captures, in particular, the mixed 
valence regime of the model, -in which electrons fluctu- 
ate from the impurity to the band and vice versa-, it 
is of utmost importance to verify that such physics is 
rigorously captured by including the Klein factors and 
showing how they take part in the transformations that 
are required for the mapping. A careful study of the 
Klein factors opens the door for future calculations of 
correlation functions using the solvable points as start- 
ing ground. The representativity of calculations that 
start from the biresonant-level model requires that such 
model flows into generic flxed points representative of the 
physics of the original model. We performed an analysis 
of the stability of deviations from the locus of solvable 
points in parameter space and established that pertur- 
bation theory on such deviations is valid. The advantage 
of starting from the mapped model is evident: closed 
analytical expressions can be derived for the full tem- 
perature crossovers of the different quantities of interest; 
this sets the approach apart from other non-perturbative 
techniques applied previously to the two-channel Ander- 
son model. 

To illustrate the versatility, power and simplicity of 
the use of the biresonant-level model to do calculations, 
we presented several results of thermodynamic quanti- 
ties that involve the impurity. Namely, we gave results 
for the impurity entropy, valence and magnetization as 
functions of temperature for different values of e and h. 
These quantities illustrate the existence of several energy 
scales in the model that signal the transitions between 
different regimes. It is remarkable that all the crossovers 
can be calculated explicitly in closed analytic expressions, 
something not possible with other types of calculations. 
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For the followup, it would be interesting to exploit 
the precise mapping (that includes all the details such 
as Klein factors, decoupled bands, etc.) to compute dy- 
namic correlation functions involving the impurity. An- 
other interesting direction would be the study of the 
finite-size crossover spectrum of the model, which would 
allow to establish instructive connections with different 
renormalization group schemes. Finally, we are in a po- 
sition to, for instance, consider the behavior of more 
than one impurity (cf. Ref. IgoI ) and, in particular, the 
case of one-dimensional two-channel Anderson lattices 
(cf. Ref. [6ll) . The latter would be interesting in order to 
determine if lattice effects are able to reverse the sign of 
the prefactor in the leading \/T term in the resistivity?^ 
which would be required in order to match the experi- 
mental results for thoriated UBeis^i^ 



(/"(r)/"(0)) ^ ^ 



(/'(r)/'(0)) - - , (71c) 



from where we read the new naive scaling dimensions: 
Wlf\ = Wss\ = W\ = [rf'] = [/I = 1/2 and [/"] = 3/2. 
For the nondiagonal Green functions, the following ones 
have long-T behaviors that do not match the expectations 
from the naive scalings: 



(v^:/t)^:^(o)) 

(^;^(r)d'(0)) 

(V':^(t)/'(o)) 

(d"(r)/"(0)) 



(^^/r)d"(0)) 
(V;;}(r)/"(0)) 
(d"(r)d'(0)) 

(rf"M/'(o)) 



1 



(72a) 
(72b) 
(72c) 
(72d) 
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APPENDIX: STUDY OF THE FINITE-FIELD 
FIXED POINTS 

The long-r asymptotics change dramatically when a 
non-zero external field is present (h ^ 0). The /' Ma- 
jorana component of / is now coupled to the rest of the 
system and we must recompute all the Green functions. 
For the diagonal ones, we find 



«}(r)<}(0)) 



1 
r 

{d"{r)d"m ^ i 



(V':/(r)V;/(0)) 



1 

T 

(d'(r)d'(O)) ^ \ 



(71a) 
(71b) 



All of them decay faster than what the naive use of the 
scaling dimensions of the fields would indicate. Thus, in 
order to obtain the RG flow of the nondiagonal couplings, 
one should, all the same, employ the naive dimensions. 
The resulting flow equations are as follows: 



dgd' 
d£ 


= -gd' , 


dgd" 
dl 


= ~gd" , 


(73a) 


dgf 

di 


= -5/' , 


dgf" 
dl 


= -35/" , 


(73b) 


dv^i 
dl 


= , 


dv^ii 
dl 


= , 


(73c) 


dA 




de 




(73d) 


H 


= , 


dl 


= , 


dr 




dh 




(73e) 


m 


= -r , 


di 


= -h . 



Now h became RG- irrelevant. However, the ratio h/T 
is exactly marginal and parametrizes the new finite-field 
Fermi-liquid fixed points. 
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